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Abstract 

Background: Magnetic Resonance Imaging (MRI) is the gold standard for detailed visualisation of spinal 
pathological and degenerative processes, but the prevailing view is that such imaging findings have little or no 
clinical relevance for low back pain. This is because these findings appear to have little association with treatment 
effects in clinical populations, and mostly a weak association with the presence of pain in the general population. 
However, almost all research into these associations is based on the examination of individual MRI findings, despite 
its being very common for multiple MRI findings to coexist. Therefore, this proof-of-concept study investigated the 
capacity of a multivariable statistical method to identify clusters of MRI findings and for those clusters to be 
grouped into pathways of vertebral degeneration. 

Methods: This study is a secondary analysis of data from 631 patients, from an outpatient spine clinic, who had 
been screened for inclusion in a randomised controlled trial. The available data created a total sample pool of 3,155 
vertebral motion segments. The mean age of the cohort was 42 years (SD 1 0.8, range 1 8-73) and 54% were women. 
MRI images were quantitatively coded by an experienced musculoskeletal research radiologist using a detailed and 
standardised research MRI evaluation protocol that has demonstrated high reproducibility. Comprehensive MRI 
findings descriptive of the disco-vertebral component of lumbar vertebrae were clustered using Latent Class Analysis. 
Two pairs of researchers, each containing an experienced MRI researcher, then independently categorised the clusters 
into hypothetical pathoanatomic pathways based on the known histological changes of discovertebral degeneration. 

Results: Twelve clusters of MRI findings were identified, described and grouped into five different hypothetical 
pathways of degeneration that appear to have face validity. 

Conclusions: This study has shown that Latent Class Analysis can be used to identify clusters of MRI findings from 
people with LBP and that those clusters can be grouped into degenerative pathways that are biologically plausible. If 
these clusters of MRI findings are reproducible in other datasets of similar patients, they may form a stable platform to 
investigate the relationship between degenerative pathways and clinically important characteristics such as pain and 
activity limitation. 
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Background 

While it is widely accepted that low back pain (LBP) is a 
bio-psycho-social condition, concern has recently been 
expressed that LBP research is moving too far from the 
biological component [1]. However, despite the large 
body of previous research attempting to identify specific 
biological reasons for LBP, there are still more questions 
than answers. The clinical consequence of this uncertainty 
is that only approximately 15% of primary care patients 
with LBP are classified with a specific pathoanatomic diag- 
nosis, with the remaining 85% being diagnosed with non- 
specific LBP [2]. 

In research, non-specific LBP is often conceptualised 
as one condition, although most clinicians and researchers 
believe that non-specific LBP is actually a heterogenous 
mix of conditions [3]. This distinction occurs because, al- 
though many methods to subgroup non-specific LBP have 
been proposed, the evidence for the validity of almost all 
these approaches remains very tentative [4]. It is also the 
case that almost all subgrouping methods that have demon- 
strated any treatment validity are not based on pathoanat- 
omy [5], the exception being the established association 
between pain centralisation and discogenic pain [6,7]. 

One way of categorising patients into diagnostic sub- 
groups is by pathoanatomy identified with Magnetic 
Resonance Imaging (MRI). MRI has increasingly re- 
placed other imaging modalities in the diagnosis of LBP 
because of the level of detail it provides. Paradoxically, 
despite the unprecedented level of detail available on 
spinal pathological or degenerative findings, the prevail- 
ing view is that they have little clinical relevance [8-10]. 
This is because these findings seem to have little or no 
association with treatment effect [11] in clinical popula- 
tions, and mostly a weak and inconsistent association 
with the presence of pain in the general population [12]. 
Furthermore, when using MRI, there is normally no way 
to distinguish age-related degeneration (believed to be 
non-painful) and 'pathological' or trauma-related degen- 
eration (believed to be painful). Similarly, it has been 
reported that normal aging of the motion segment and 
pathological degeneration cannot be distinguished at a 
histological level [13], suggesting that the biological 
process is the same in normal aging but is accelerated in 
the pathological state. Therefore, the biological age of a 
vertebral motion segment is not necessarily the same as 
the actual chronological age of the person. 

However, investigation into associations between MRI 
findings and pain are complicated by numerous MRI find- 
ings being present at the same time. For example, verte- 
bral endplate signal changes and vertebral disc herniations 
almost always co-exist with other degenerative disc find- 
ings, such as reduction of height and reduction of the signal 
intensity of the disc [14-17]. Despite this, most previous re- 
search has focused solely on the association between single 



imaging findings and pain or other clinical outcomes. It is 
only recendy that researchers have started to explore the 
potential for answers within this complexity, for example 
Cheung and co-workers[18], who have reported a positive 
correlation between the sum of degenerative disc MRI find- 
ings and low back pain. 

Classification systems exist for categorising the de- 
tailed information on the physiological condition of a 
motion segment that MRI provides. For example, the 
Pfirrmann classification [19] is a descriptive tool for clas- 
sifying intervertebral disc degeneration. However, such 
classification systems also typically focus on findings re- 
lated to one or two anatomic entities, such as the height 
and signal intensity of the intervertebral disc, leaving 
un-described other relevant entities such as vertebral 
endplate irregularity, vertebral endplate signal changes 
and disc herniations. 

It would be useful if there were methods to better 
model the multivariable relationships between clusters 
of MRI findings, as these might provide a clearer under- 
standing of how pathoanatomical processes are expressed 
across all the structures associated with a vertebral seg- 
ment. Furthermore, the identification of such processes, 
and an ability of such clusters to stage the biological age 
of a vertebral segment, might allow better insight into the 
relationship between degenerative processes and clinical 
characteristics such as pain and activity limitation. 

Therefore, the aims of this study were to investigate: 
(i) how multiple MRI findings from lumbar spine motion 
segments cluster together in patients with chronic LBP, 
and (ii) to classify these clusters of MRI findings into 
hypothetical pathoanatomical pathways. 

Methods 

This was a 'proof of concept' study designed to test a novel 
analytical method and set hypotheses about pathoanat- 
omical pathways of vertebral segment degeneration. It 
used cross-sectional data and a probabilistic (Bayesian) La- 
tent Class Analysis approach to the multivariable cluster- 
ing of MRI findings. 

Study sample 

This study is a secondary analysis of data from a cohort 
of patients who were potential participants in a ran- 
domised controlled trial [20]. All participants attended the 
same publically-funded outpatient spine clinic (The Spine 
Centre of Southern Denmark) where they had been re- 
ferred from the primary care sector for a multidisciplinary 
evaluation. In this clinical setting, from June 2006 to June 
2008, MRI was routinely performed on all patients (who 
had no contraindications for MRI) who met the following 
criteria: (a) LBP or leg pain of at least 3 on an 11 -point 
Numerical Rating Scale, (b) duration of current symptoms 
from 2 to 12 months, and (c) age above 18 years. 
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Data collection 

MRI was performed with a 0.2 T MRI-system (Magnetom 
Open Viva; Siemens AG, Erlangen, Germany). A body 
spine surface coil was used for imaging of the lumbar re- 
gion, with the study subjects in the supine position. The 
imaging protocol consisted of the following sequences: 

• Localizer sequence, 40/10/40 (TR/TE/flip angle), two 
coronal and three sagittal images in orthogonal planes 

• Sagittal Tl-weighted spin echo, 621/26 (TR/TE), 
144 x 256 matrix, 300 mm FOV, and 11 slices 4 mm 
wide, distance factor 0.20. 

• Sagittal T2-weighted turbo spin echo, 4609/134 
(TR/effective TE), 210 x 256 matrix, 300 mm FOV, 
and 11 slices 4 mm wide, distance factor 0.20. 

• Axial Tl-weighted spin echo, 720/26 (TR/TE), 192 x 
256 matrix, 240 mm FOV, and 15 slices 5 mm wide, 
distance factor 0.25. 

• Axial T2-weighted turbo spin echo, 6415/134 (TR/ 
effective TE), 180 x 256 matrix, 250 mm FOV, and 
15 slices 5 mm wide, distance factor 0.25. 

• Axial images were performed on the three lower 
lumbar levels. If herniations were present at higher 
lumbar levels, relevant supplementing axial series 
were performed. 

The images were quantitatively coded using a detailed 
and standardised research MRI evaluation protocol [21,22]. 
The evaluation was performed by an experienced musculo- 
skeletal research radiologist who was blinded to any par- 
ticipant information other than name, age and sex. 

Variables of interest 

The following MRI variables were included in the current 
study: intervertebral disc height, disc signal intensity, disc 
herniations, size and type of vertebral endplate signal 
changes (VESC), irregularity of the vertebral endplate, 
osteophytes and spondylolisthesis. Each MRI variable was 
coded on two to seven different pathoanatomical stages 
(categories), as detailed in Additional file 1. 

Previous testing of this MRI evaluation protocol when 
used by the same radiologist has shown substantial to al- 
most perfect inter- and intra-observer reliability for the 
VESC variables (Kappa ranging from 0.73 to 1.0), and 
moderate to substantial reliability for other endplate 
related findings (Kappa ranging from 0.52 to 0.72) [21]. 
Evaluation of intervertebral disc-related changes has 
shown moderate to almost perfect reliability (Kappa ran- 
ging from 0.59 to 0.97) [22]. Anterolisthesis was evalu- 
ated according to Meyerding classification system [23]. 

Statistical analyses 

MRI findings were clustered using the multivariable stat- 
istical program 'SNOB' (Monash Data Mining Centre, 



Monash University, Melbourne, Victoria, Australia), which 
is a form of Latent Class Analysis. SNOB uses probabilis- 
tic mixture modelling and the minimum message length 
principle to determine the optimum class structure (clus- 
ters) within the data. Minimum message length is an auto- 
mated statistical method to determine which of the 
possible models of class structure and explained variance 
are the most parsimonious and explanatory. Latent Class 
Analysis has a number of advantages over traditional clus- 
ter analysis techniques, including: greater classification 
accuracy, the ability to manage variables of all data types 
(dichotomous, ordinal and continuous), and a tolerance of 
missing data [24-26] . 

Analysis occurred at the level of individual vertebral 
segments, as lumbar segments within an individual person 
may each have a different 'pathoanatomical age'. There- 
fore, every participating person contributed five lumbar 
vertebral segments to the analysis. 

The Latent Class Analysis identified the optimal num- 
ber of clusters and also the cluster membership for each 
vertebral segment. We used this information to deter- 
mine the proportion of vertebral segments within each 
cluster that displayed each coding category on each MRI 
variable and the proportion of vertebral levels (Ll/2, L2/ 
3, L3/4 L4/5 and L5/S1) within each cluster. These data 
were calculated and graphed using Excel 2008 for Mac ver- 
sion 12.2.8 (Microsoft Corporation, Redmond, WA, USA). 

Then two pairs of researchers, each containing an ex- 
perienced MRI researcher, independently categorised the 
clusters into hypothetical pathoanatomic pathways based 
on the face validity of known histological changes of 
discovertebral degeneration. These histological changes 
had been identified by an (unpublished) electronic litera- 
ture search (PubMed and Medline) and review. The 
pathways categorised by the pairs of researchers were 
compared and any differences resolved by discussion. 

A post-hoc calculation of the mean (SD) chronological 
age of the motion segments within each cluster was also 
performed. This metric was used to test if the chrono- 
logical age of each cluster challenged or supported the 
concept of 'pathoanatomical age' inherent in the pathways. 

Ethics 

This analysis was based on existing data collected for a 
randomised controlled trial [20] approved by the Ethics 
Committee for the Region of Southern Denmark (ap- 
proval # S-VF-20060111), registered in ClinicalTrials.gov 
(Identifier # NCT00454792) and performed following 
the Declaration of Helsinki principles. 

Results 

Population 

MRI findings and basic demographic data were available 
on 631 patients who had been screened for inclusion in 
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the original studies. All were included in the current 
study cohort, which created a total sample pool of 3,155 
vertebral motion segments. The mean age of the cohort 
was 42 years (SD 10.8, range 18-73) and 54% were 
women. Detailed descriptive data of the MRI findings at 
a whole group level are tabulated in Additional file 1. 

Clusters of MRI findings 

The Latent Class Analysis identified twelve clusters of 
MRI findings. One cluster, characterised by no abnormal 
MRI findings, contained 52% of the 3,155 motion seg- 
ments and represented the normal, pre-degenerative 
state. The second largest cluster contained almost 15% 
of the motion segments and was characterised by mo- 
tion segments with reduced disc height, reduced disc 
signal intensity, disc bulges and a minor degree of disc 
protrusions and high intensity zones. Sixty percent of 
these changes were located at the two lowest spinal levels. 
The rest of the clusters were smaller, with each containing 
less than 8% of the motion segments. Each cluster and its 
distribution of MRI findings is shown diagrammatically in 
Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, and 12. More detailed 
descriptive data of the MRI findings at a cluster level are 
tabulated in Additional file 2. 



Hypothetical pathoanatomic pathways 

The hypothetical pathoanatomic pathways derived from 
the content analysis of the clusters were: (i) two clusters 
representing progressive stages of disc degeneration in the 
lower lumbar levels; (ii) four clusters representing progres- 
sive stages of disc herniations and VESC in the lower lum- 
bar levels; (iii) two clusters containing progressive endplate 
changes and disc degeneration at the upper lumbar levels 
only; (iv) two cluster with VESC only, one cluster each for 
changes in the upper and the lower endplates; and lastly, 
(v) one cluster containing osteophytes at the upper lumbar 
motion segments. The pathoanatomic pathways and their 
clusters are illustrated in Figure 13. 

Once the pathoanatomic pathways had been catego- 
rised, the mean age of the motion segments in each clus- 
ter was added to Figure 13 as a post-hoc analysis. The 
average age of the vertebral segments within each cluster 
supported the notion of a developmental model, as the 
chronological age in each cluster consistently increased 
across each pathoanatomic pathway as degeneration de- 
veloped. Although the chronological age and biological 
age of vertebral segments may not be the same due to 
pathological degeneration, the observed age progression 
across the biological pathways adds credibility to the 
construct validity of this classification. 



Cluster 1: 'No findings' 



90% 
80% 
70% 
60% 
50% 
40% 
30% 
20% 
10% 
0% 



Disc findings 



Vertebral findings 



L4 
Il5 



/ # i i J i & f & f f f f f & & J 
////// / / / ✓ / / / / / / /. 

f // / ////////// * 4 



V 



1 Gat. 1 1 Cat. 2 □ Cat. 3 □ Cat. 4 □ Cat. 5 □ Cat. 6 

Figure 1 Prevalence of MRI findings within Cluster 1. This cluster contained 51 .6% of the 3155 vertebral motion-segments, (a) The vertical 
bars on the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars 
on the graph do not show the 'category 0' as that represents the presence of only normal findings. For pathologies with only one category it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC the categories refer to increasing severity or 
size of pathology with 'Category V being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 
2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC, the categories refer to the following types: 1 =Type I; 2 =Type 
II; 3 =Type III; 4= mixed Type l/ll; 5 = mixed Type ll/lll; 6 = mixed Type \/\\\. For 'anterolisthesis', the categories refer to the Meyerdings 
classification. Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the 
motion segments was 39 years (SD 10) and the proportion of women was 58%. (c) VESC = vertebral endplate signal change; 
Cat. = coding categories; EP = endplate. 
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Cluster 2: 'Moderate disc degeneration' 
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Figure 2 Proportion of categories of MRI findings for Cluster 2 containing 14.9% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the 'Category 0', as that represents the presence of only normal findings. For pathologies with only one category, it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC' the categories refer to increasing severity or 
size of pathology with 'Category V being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC, the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/ll; 5 = mixed Type ll/lll; 6 = mixed Type l/lll. For 'anterolisthesis' the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
46 years (SD 1 1) and the proportion of women was 55%. (c) VESC = vertebral endplate signal change; Cat. = categories; EP = endplate. 



Cluster 3: 'Severe disc degeneration' 
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Figure 3 Proportion of categories of MRI findings for Cluster 3 containing 7.7% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the 'Category 0', as that represents the presence of only normal findings. For pathologies with only one category, it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC, the categories refer to increasing severity or 
size of pathology with 'Category V being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC, the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/ll; 5 = mixed Type ll/lll; 6 = mixed Type l/lll. For 'anterolisthesis' the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
47 years (SD 10) and the proportion of women was 47%. (c) VESC = vertebra! endplate signal change; Cat. = categories; EP = endplate. 
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Cluster 4: 'Disc degeneration and herniations' 
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Figure 4 Proportion of categories of MRI findings for Cluster 4 containing 5.5% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the 'Category 0', as that represents the presence of only normal findings. For pathologies with only one category, it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC, the categories refer to increasing severity or 
size of pathology with 'Category 1' being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/l I; 5 = mixed Type ll/lll; 6 = mixed Type l/l 1 1. For 'anterolisthesis', the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
41 years (SD 10) and the proportion of women was 48%. (c) VESC = vertebral endplate signal change; Cat. = categories; EP = endplate. 



Cluster 5: 'Disc degeneration, herniations, osteophytes and large VESC 
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Figure 5 Proportion of categories of MRI findings for Cluster 5 containing 5.4% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the 'Category 0', as that represents the presence of only normal findings. For pathologies with only one category, it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC, the categories refer to increasing severity or 
size of pathology with 'Category V being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC, the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/l I; 5 = mixed Type ll/lll; 6 = mixed Type l/lll. For 'anterolisthesis' the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
48 years (SD 9) and the proportion of women was 51%. (c) VESC = vertebral endplate signal change; Cat. = categories; EP = endplate. 
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Cluster 6: 'Disc degeneration, herniations, osteophytes and small VESC 
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Figure 6 Proportion of categories of MRI findings for Cluster 6 containing 3.3% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the 'Category 0' as that represents the presence of only normal findings. For pathologies with only one category, it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC, the categories refer to increasing severity or 
size of pathology with 'Category 1' being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC, the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/l I; 5 = mixed Type ll/lll; 6 = mixed Type l/l 1 1. For 'anterolisthesis', the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
47 years (SD 9) and the proportion of women was 53%. (c) VESC = vertebral endplate signal change; Cat. = categories; EP = endplate. 



Cluster 7: 'Osteophytes present at upper levels' 
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Figure 7 Proportion of categories of MRI findings for Cluster 7 containing 3.0% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the Category 0', as that represents the presence of only normal findings. For pathologies with only one category, it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC, the categories refer to increasing severity or 
size of pathology with Category V being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC, the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/l I; 5 = mixed Type ll/lll; 6 = mixed Type l/lll. For 'anterolisthesis', the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
45 years (SD 9) and the proportion of women was 64%. (c) VESC = vertebral endplate signal change; Cat. = categories; EP = endplate. 
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Figure 8 Proportion of categories of MRI findings for Cluster 8 containing 2.1% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the 'Category 0' as that represents the presence of only normal findings. For pathologies with only one category, it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC, the categories refer to increasing severity or 
size of pathology with 'Category V being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC, the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/l I; 5 = mixed Type ll/lll; 6 = mixed Type l/lll. For 'anterolisthesis', the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
35 years (SD 9) and the proportion of women was 28%. (c) VESC = vertebral endplate signal change; Cat. = categories; EP = endplate. 



Cluster 9: 'VESC - present at upper endplates' 
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Figure 9 Proportion of categories of MRI findings for Cluster 9 containing 2.0% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the 'Category 0', as that represents the presence of only normal findings. For pathologies with only one category it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC, the categories refer to increasing severity or 
size of pathology with 'Category V being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC, the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/l I; 5 = mixed Type ll/lll; 6 = mixed Type l/lll. For 'anterolisthesis', the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
47 years (SD 7) and the proportion of women was 34%. (c) VESC = vertebral endplate signal change; Cat. = categories; EP = endplate. 
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Cluster 10: 'Disc degeneration, herniations and osteophytes' 
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Figure 10 Proportion of categories of MRI findings for Cluster 10 containing 1.7% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the 'Category 0', as that represents the presence of only normal findings. For pathologies with only one category it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC, the categories refer to increasing severity or 
size of pathology with 'Category V being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC, the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/l I; 5 = mixed Type ll/lll; 6 = mixed Type l/lll. For 'anterolisthesis', the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
43 years (SD 10) and the proportion of women was 50%. (c) VESC = vertebral endplate signal change; Cat. = categories; EP = endplate. 
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Cluster 1 1 : 'VESC - present at lower endplates' 
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Figure 1 1 Proportion of categories of MRI findings for Cluster 1 1 containing 1.6% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the 'Category 0', as that represents the presence of only normal findings. For pathologies with only one category, it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC, the categories refer to increasing severity or 
size of pathology with Category V being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC, the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/l I; 5 = mixed Type ll/lll; 6 = mixed Type l/lll. For 'anterolisthesis', the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
44 years (SD 1 0) and the proportion of women was 42%. (c) VESC = vertebral endplate signal change; Cat. = categories; EP = endplate. 
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Figure 12 Proportion of categories of MRI findings for Cluster 12 containing 1.1% of the 3155 motion-segments, (a) The vertical bars on 
the graph represent the proportion of vertebral motion segments present within each coding category for each MRI finding. The bars on the 
graph do not show the 'Category 0', as that represents the presence of only normal findings. For pathologies with only one category, it 
represents the presence of the pathology. For 'disc signal intensity', 'disc height' and 'size of VESC, the categories refer to increasing severity or 
size of pathology with 'Category V being the least severe. For 'herniations', the categories represent the following types: 1 = focal protrusion; 

2 = broad-based protrusion; 3 = extrusion; 4 = sequestration. For 'type of VESC the categories refer to the following types: 1 =Type I; 2 = Type II; 

3 =Type III; 4 = mixed Type l/l I; 5 = mixed Type ll/lll; 6 = mixed Type l/l 1 1. For 'anterolisthesis' the categories refer to the Meyerdings classification. 
Detailed information on proportions of categories of MRI findings is described in Additional file 2. (b) The mean age of the motion segments was 
37 years (SD 12) and the proportion of women was 41%. (c) VESC = vertebral endplate signal change; Cat. = categories; EP = endplate. 



Discussion 

In this study, comprehensive MRI findings of lumbar 
vertebral motion segments were grouped using Latent Class 
Analysis into twelve clusters. Using content analysis, the 
identified clusters were then qualitatively grouped into five 
different biological pathways of degeneration that appear to 
have face validity. 

This data-driven approach of the analysis of MRI find- 
ings is novel in LBP and our proof-of-concept study 
shows that this methodological approach produces plau- 
sible results suitable for further investigation. Two previ- 
ous studies on LBP have used a multivariable approach to 
exploring the relationship between MRI findings and clin- 
ical variables, but in contrast to the current study, their 
multivariable analysis clustered the clinical variables and 
related these to single MRI findings. Takatalo et al. [27,28] 
used Latent Class Analyses to group multiple clinical vari- 
ables (including pain, activity limitation and care seeking) 
from 468 people into five clusters which were then investi- 
gated for their association with individual MRI findings. It 
is not yet known whether the most revealing investigations 
of the association between MRI findings and clinical char- 
acteristics will involve the clustering of findings from both 
domains, only one domain or neither domain. In other 
areas of research, such as neurology and biological psych- 
iatry, Latent Class Analysis has been used with MRI 



findings to investigate the validity of diagnostic subgroups 
or the accuracy of diagnostic tests [29-32]. 

Although the emphasis in the current study was 
on method, the identified MRI clusters appear to be bio- 
logically plausible and could be grouped into pathways 
of degeneration. These suggested pathways were based 
on the available evidence about the histological progres- 
sion of disco-vertebral degeneration for different age 
groups [13] and on classification systems of degeneration 
based on histological [33] and MRI findings [19]. The 
chronology of the degenerative pathways was also in- 
formed by knowledge of the potential reversibility of 
some MRI findings (for example, an intervertebral disc 
herniation is potentially reversible but reduced disc height 
is not) and also, if the distribution of findings were mainly 
located in the lower or the upper lumbar motion segments. 
However, as there are few studies describing the precise as- 
sociation between the histology of degeneration and spe- 
cific MRI findings, the evidence on which to base such 
grouping of MRI clusters is incomplete. Therefore, the 
clusters and degenerative pathways should be interpreted 
with caution and viewed as hypothesis-setting only. 

MRIs of the lumbar spine can provide researchers 
with highly detailed information on spinal pathoanatomy, 
especially when quantified via comprehensive research 
protocols [21,22]. However, the volume of this detail can 
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Figure 13 Degenerative pathways of the disco-vertebral complex. 'Age' is the average chronological age of the vertebral segments in each 
cluster, 'SD' is the standard deviation and 'VESC is vertebral endplate signal change. The clusters were divided into the following five pathways: 
(i) Progressive stages of disc degeneration in the lower lumbar levels, (ii) Progressive stages of disc herniations and VESC in the lower lumbar 
levels, (iii) Progressive endplate changes and disc degeneration at the upper lumbar levels only, (iv) VESC only at either the upper or the lower 
endplates. (v) Osteophytes at the upper lumbar motion segments. 



present a challenge for data analysis. For example, in our 
study, for each disco-vertebral level, there were more than 
50 variables, each with up to seven response options. 
Traditional approaches to data reduction, such as regres- 
sion analysis, can model the dominant pattern of associ- 
ation between MRI variables but do not accommodate the 
presence of multiple patterns that may reflect different 
stages of degeneration or different pathological pathways. 
Clustering techniques, such as Latent Class Analysis, have 
the advantage of allowing for that complexity. 

We interpret these results as indicating that this re- 
search method should be further tested in different 
datasets to investigate if similar clusters are consistently 
replicated or if such clusters are highly sample-specific. 
If these clusters are reproducible, they could provide a 
platform from which to investigate associations between 
multivariate MRI findings and clinical findings, such as 
pain, activity limitations or other clinically relevant out- 
comes. Similarly, a further extension of this work would 
be to include MRI findings of structures other than the 
disco-vertebral complex, such as MRI findings of the 



vertebral canal, zygapophyseal joints and para-spinal soft 
tissue structures. 

The strengths of this study are that a new statistical ap- 
proach was applied to a large sample of MRI findings that 
had been rigorously described by an experienced research 
radiologist using a data extraction protocol with high re- 
producibility, and the results were cautiously interpreted. 
A weakness of the study is that its methodological focus 
means that the content of the identified clusters and bio- 
logical pathways are preliminary and require extensive fur- 
ther study before their clinical validity, if any, can be 
determined. 

Conclusions 

This study has shown that Latent Class Analysis can be 
used to identify clusters of MRI findings from people with 
LBP and that those clusters can be grouped into biological 
pathways that have some face validity. If, when applied to 
other datasets of similar patients, this research method 
identifies reproducible clusters of MRI findings, these may 
form a stable platform to investigate the relationship 
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between degenerative pathways and important clinical 
characteristics such as pain and activity limitation. 

Additional files 



Additional file 1: MRI findings and their distribution in the whole 
sample of vertebral segments. Detailed descriptive data of the MRI 
findings at a whole group level are tabulated in Additional file 1. 

Additional file 2: The distribution (%) of the MRI findings in the 12 
clusters. Detailed descriptive data of the MRI findings at a cluster level 
are tabulated in Additional file 2. 
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